Polar magneto-optical Kerr effect for low-symmetric ferromagnets 
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The polar magneto-optical Kerr effect (MOKE) for low-symmetric ferromagnetic crystals is investigated 
theoretically based on first-principle calculations of optical conductivities and a transfer matrix approach for 
the electrodynamics part of the problem. Exact average magneto-optical properties of polycrystals are de- 
scribed, taking into account realistic models for the distribution of domain orientations. It is shown that for 
low-symmetric ferromagnetic single crystals the MOKE is determined by an interplay of crystallographic bire- 
fringence and magnetic effects. Calculations for single and bi-crystal of hep { 1 120) Co and for a polycrystal of 
Cr02 are performed, with results being in good agreement with experimental data. 

PACS numbers: 78.20.Bh,78.20.Ci 



5h 

43 



o 
o 



> 

in 
o 
in 
o 
-i— > 



T3 

G 
O 

o 



X 



I. INTRODUCTION 

The magneto-optical Kerr effect (MOKE) is a versatile 
method to probe magnetic properties of thin films. Advanced 
by the rapid developments in crystallographic growth tech- 
niques, a variety of low-symmetric crystalline surfaces have 
been subject to MOKE measurements in the last decades. 
This has led to systematic investigations of magneto-optical 
anisotropy effects^. 

State of the art theoretical approaches to investigate 
the MOKE are based on first-principle calculations of di- 
electric tensors in the framework of the Kubo-Greenwood 
formalism 2 - 3 as suggested by Wang and Callaway 4 . The 
MOKE is obtained from a dielectric tensor by means of an 
approximative analytic expression 
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derived originally by Argyres in 1955 5 . \|/ denotes the Kerr 
rotation and % denotes the Kerr ellipticity. 

This approach requires in general that the dielectric tensor 
has symmetry 



(2) 



There have been theoretical attempts to extend the approach 
to low-symmetric systems, however so far the complete elec- 
trodynamics calculation for low symmetric dielectric tensors 
has not been considered. 

There are many interesting ferromagnets that have a low 
symmetry, e.g. Cr02, hep (1120) Co and FePt grown in the 
(010) direction. All of these systems have two different crys- 
tallographic axis in the surface plane, so beside their magneto- 
optical activity they exhibit crystallographic birefringence. 

In this paper we show that for such crystals it is important 
to consider the complete optical response including birefrin- 
gence and magnetic effect in order to describe correctly the 
polar MOKE. Further, we show that the optical response is 
qualitatively different for single- and polycrystals and finally, 



for polycrystals it sensitively depends on the ordering of crys- 
tallographic domains. We calculate the MOKE of hep (1 120) 
Co and of (010) Cr02- For Co we show that the previous in- 
terpretation of experimental data of anisotropic polar MOKE 1 
in terms of a manifestation of magneto-crystalline anisotropy 
remains valid. 

The paper is organised as follows. In the subsequent section 
we describe our approach to the complete calculation of the 
electrodynamics problem by means of transfer matrix meth- 
ods. Theoretical description of ellipsometry measurements 
for single- and polycrystals is given in Sec. [ni| In Sec. IIVI 
we discuss first-principle calculations of optical conductivi- 
ties. Space-time symmetry of Co and Cr02 crystals is de- 
scribed in Sec.|V] The calculated optical response of Co and 
QO2 is presented in Sec. IV II and Sec. IV 111 respectively. In 
Sec. I Villi a summary and conclusions are given. 



II. TRANSFER MATRIX METHODS 

The optical response of a finite system of layers to an 
incident plane wave can be described by transfer matrix 
methods^i 7 ^. The description is valid if the magnetic perme- 
ability is unity and the wavelength of the light is large com- 
pared to the microscopic structure of materials and also large 
compared to interface roughness. In the most general case a 
system with n boundaries is described by a regular set of An 
linear equations that determines the complex amplitude vec- 
tors of all plane waves in all media. We briefly describe the 
method. 

We first choose a coordinate system such that the z-axis is 
the surface normal and the scattering plane is spanned by the 
z-axis and the y-axis. In the half space of the incident and 
reflected wave Fresnel's secular equation reads 



7 CO 

I^E+^rEE = 0. 



(3) 



We substitute r i— » — r and define q := j-k v and k := j-k 7 . This 

C " CO > CO * 



2 



gives 



, k 4 are the roots of the fourth order polynomial in k 



q = Vesin-r) 

k = ±\/ecos-r} =: ±ko, 



(4) 



where $ is the incident angle. This gives an ansatz for the 
wave 
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where E' n is the known amplitude vector of the incident wave 
and the complex amplitude vector of the reflected wave satis- 
fies 



rprefl _ 1 Tprefl 
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leaving two free parameters E x and E r y e ^ . For other media 
the most general plane wave solution to Maxwell's equations 
is a combination of four independent waves. In the case of a 
scalar medium it is 



where 
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and the x- and y-components of E l and E 2 are independent. 
In the case of a tensor medium it is 

E = a l n l e^ y + k[z - at '> + ...+ fl Ve ! '^ + * Vo *) (9) 

with four free parameters a 1 , . . . ,a 4 satisfying 

(10) 



Det 



e« - q 



Eyy — k 2 e yz + qk 
e zy + qk e zz — q 2 



(11) 



and the vectors n 1 , . . . ,n 4 are associated kernels. 

In the half space on the backside of the layers two waves 
can always be discarded. For transparent medium these are 
two backward travelling waves, for an absorbing medium 
these are two exponentially decaying waves. 

In our case (a bulk metallic system with no intermediate 
layer) we have only an absorbing tensor half space and the 
ansatz for the waves in the responding system reduces to 



E = flWfo** 1 *-™) +a 2 n 2 e i ^ y+k2z - (0t \ 



(12) 



where k l and k 2 are the roots that have negative imaginary 
parts (negative z-direction corresponds to forward travelling 
waves). 

Stressing the assumption of unity magnetic permeability, 
four independent boundary conditions follow from Maxwell's 
equations stating that 



E X) E y , d z E x and iqE z — d z E y 



(13) 



are continuous. 

Substituting the ansatz, Eq. (|5} and Eq. (I12> in the boundary 
conditions, we get 
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This is a regular system of four linear equations. Stressing 
Eq. l|6) and Eq. d 1 01 its solution determines the complex am- 
plitudes vectors of all waves. 



We have written a numerical implementation of the most 
general case of a transfer matrix approach (based on standard 
LAPACK 9 routines and polynomial solver 10 ). It is described 
in detail in Ref.fTTl 



III. ELLIPSOMETRY FOR SINGLE- AND 
POLYCRYSTALS 

The state of polarization of a plane wave is conveniently 
described by Stokes parameters^ 



S = 



( E X Ex + EyEy \ 
ii'rfir E yEy 
ExEy ~f~ E \Ey 

\i(E x E y -E x E y )J 



(15) 



where E = (E x ,E y ) is the complex amplitude vector of the 
plane wave in the coordinate system of the polarization state 
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analysis. 

The state of polarization of a set of incoherent plane waves 
that add by their intensities is described by the sum of their 
Stokes parameters. Both for a single wave and for an incoher- 
ent wave, the rotation angle of the polarization ellipse \\t and 
its ellipticity % are related to the Stokes parameters by 

5? 

tan2\|/= -f (16) 
Si 

and 

sin2x=^===. (17) 

In the general case the polarization ellipse is the intensity be- 
hind an analyser for all positions. Only in the special case of 
a single wave is this equivalent to the curve that is drawn by 
the tip of the electric field vector. 

The optical response of a polycrystal can be described by 
the sum over Stokes parameters of single crystalline domains 
weighted by surface areas of the domains and intensities shin- 
ing on them 12 . The sum extends over all domains that are 
illuminated in the experiment. The approach is valid if single 
crystalline domains are large compared to the wavelength. We 
can calculate Stokes parameters for polycrystals by summing 
over Stokes parameters obtained from transfer matrix calcula- 
tions for single crystals. 

IV. FIRST-PRINCIPLE CALCULATIONS OF OPTICAL 
CONDUCTIVITIES 

We briefly describe the calculation of optical constants by 
means of first-principle calculations. Our approach is basi- 



cally standard unless we evaluate the Kubo-Greenwood for- 
mula directly without Kramers-Kronig transformation and 
analytical continuation (see also Ref. ll3l) . 

In this section we consider the optical conductivity tensor o 
rather than the corresponding dielectric tensor e. The quanti- 
ties are related by the identity 

e a p(ffl) = 5 p + i§a a p(co). (18) 

In general, intra-band, as well as direct and indirect inter- 
band transitions, contribute to the optical conductivity. Spins 
may flip (for magnetic dipole transitions) or stay constant (for 
electric dipole transitions ) during excitations. It is a com- 
mon practice to account only for the contribution of electric 
dipole (non-spin-flip) direct inter-band transitions by means 
of ab initio methods while treating the contribution of intra- 
band transitions by a phenomenological Drude term 

1 + CO z T z 

and neglecting all other contributions 4 *! 4 -!!!^. A broad vari- 
ety of linear optical and magneto-optical effects in metals as 
well as in semiconductors have been successfully d escribe d in 
the framework of this approximation, see e.g. Refs. ll7ll8l and 
references therein. In the transition metals, intra-band transi- 
tions turn out to be important in the range from eV up to 0.5 
e\4i^. It is shown in Ref.0that a corresponding Drude contri- 
bution is negligible for energies larger than 1 e V in the case of 
Ni. Throughout this work we neglect any phenomenological 
Drude contribution. 

The Kubo-Greenwood expression for the contribution of 
direct inter band transitions to the optical conductivity reads 14 
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where the indices / and n denote the spin and all band quan- 
tum numbers for the occupied and empty states respectively 
and k is the quasi momentum running through the Brillouin 
zone, Ef is the Fermi energy. The symbol II" (k) , a = x,y,z 
denotes the matrix elements of the momentum operator given 
below by Eq. i22\ . and CO,,/ (k) is the energy difference between 
the involved states, 



n 



(21) 



Finally, t(co) is a phenomenological relaxation time. 
Throughout this work we use a constant relaxation time of 



0.136 eV. The results of this paper are insensitive to the actual 
choice of this value. 

Together with the energy differences CO,,/ (it), the matrix ele- 
ments of the momentum operator are obtained from the under- 
lying band structure calculation by evaluating the expression 



U,„(k)= d'r Mf\(k,r) 



p + — -[oxVl/(r)] 

4mc z 



\\r„(k,r) 



(22) 

Here \\i n (k,r) is the Bloch wave function with quantum num- 
bers as described above, p — —ihV and V(r) is a crystal po- 
tential. State-of-art works on ab initio calculated optical con- 
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stants neglect the spin-orbit term in the expression for the ma- 
trix elements of the momentum operator, Eq. J22i . This has 
been found to be a good approximation, see, e.g., Ref. 4. We 
follow this approach. 

Expression J20b may be computed directly or via sym- 
metrized limit expressions requiring Rramers-Kronig trans- 
formations and analytical continuation to finite relaxation 
times. We recently discussed advantages and disadvantages 
of both approaches that become important when the conduc- 
tivity tensor has low symmetry 1 3 . In the present paper Expres- 
sion j20t is computed directly. 

For electronic structure calculation we use a relativistic full- 
potential linear muffin-tin orbital (FP-LMTO) code. The code 
is described in detail in Ref. 20. A discussion of the treatment 
of spin-orbit coupling by means of the second variational step 
can be found in Ref.l2ll 

V. SYMMETRY CONSIDERATIONS 



Next we consider the symmetry properties of the same crys- 
tals but without magnetism. Co has the well known point 
group 6/mmm and irreducible form of the dielectric tensor 
without magnetism is 

/e» 0\ 

E = Eyy . (28) 

\0 Zyy) 

The &O2 crystal without magnetism is non-symmomorphic. 
It has space group P42/mmm. Evaluation of Neumann's prin- 
ciple is standard for pure point group operators. For symmetry 
operators 6 = aoT that are a combination of a point group op- 
erator a and the translation operator T (the 4-fold screw axis 
4xoT(c/2,0,0) in our case) Neumann's principle can be eval- 
uated by stressing the invariance of the dielectric tensor under 
arbitrary translations. 

ro£or'=£. (29) 



We have used standard space-time symmetry analysis— to 
find the irreducible forms of the dielectric tensors of hep Co 
with magnetisation along (1120) and of CrC<2 with magneti- 
sation along (010). The crystal structures are shown in Fig.^ 
We find space-time point groups mpil and 2/m for Co and 
C1O2 respectively. Making the coordinate systems explicit, 
irreducible sets of point group operations can be chosen as 
identity, 2-fold rotation around z followed by space inversion 
and 2-fold rotation around y followed by time inversion for Co 
and identity, space inversion and 2-fold rotation around x fol- 
lowed by time inversion for CrC<2. Standard symbols are 1, 2z, 
2y and 1, 1, 2x respectively. 

Irreducible space-time symmetries of the dielectric tensors 
follow by Neumann's principle which states that 



ooeoo 



(23) 



has to be satisfied for any symmetry operator a. For classical 
point group operators the respective matrix equation can be 
evaluated. For non-classical operators a = sox composed of 
a classical operator s and the time inversion operator T, Eg. 1231 
can be brought in matrix form by stressing the equivalence of 
time inversion and magnetisation reversal, 

To£(M)or' =e(-M) (24) 

and Onsager's relation, 

e(-Af) = e r (M), (25) 

where T denotes the transpose. 
For the Co crystal we find 



For CrCh we have 
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(26) 



(27) 



We find the irreducible form of the dielectric tensor without 
magnetism is the same as for Co. 

Next we consider the expansion of the dielectric tensor in 
powers of the magnetisation and stress the following symme- 
try properties: The zero order contribution has symmetry of 
the non-magnetic crystal. Magnetic contributions of odd or- 
der have space-time symmetry of the magnetic crystal and are 
anti symmetric. Magnetic contributions of even order have 
space-time symmetry of the magnetic crystal and are sym- 
metric. Anti symmetry respectively symmetry property of odd 
and even order magnetic contributions are arrived at in general 
by applying Onsager's relation to the expansion. 

We find that up to second order in the magnetisation the 
expansion has the symmetry, for Co, 
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and for CrC»2 
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(31) 

Results of standard electronic structure calculations are 
for both systems tensors of the formi2i22i2i (see also 
Secs. lVIAlVTTAl 



(32) 



This has an important implication. It means that second order 
magnetic contribution (which would appear as e.g. a differ- 
ence between e yy and £,,) is either absent in both systems or 
not resolvable with standard electronic structure calculation. 
There is no reason why second order magnetic contribution 
should be absent. So basically the conclusion is that it is not 
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FIG. 1: Crystal structures of a) (1120) ferromagnetic hep Co and 
b) (010) ferromagnetic Cr02- Coordinate systems are as used in 
magneto-optics calculation. 



resolvable with standard electronic structure calculations. We 
discuss this in more detail in Sec. lVIl 

For the case of CrC>2 we conclude that z xz is actually zero 
in the first order magnetic contribution, however it might still 
be present in the third order. 



VI. POLAR MOKE OF { 1 120) HCP CO 
A. Optical conductivity 
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FIG. 2: Calculated optical conductivity tensor of hep Co with mag- 
netisation direction (1120). Quantities are shown in a range where 
only direct inter band transition are important. 



rameters of our calculation (e. g. basis set, k-point mesh) the 
variation of a yy and a,, is typically larger. So we have to con- 
clude within the error of our calculation a yy and a zz are equal. 
The conclusion is that second order magnetic contribution can 
not be resolved with standard electronic structure calculation. 



We have calculated the optical conductivity tensor of hep 
(1 120) Co. A hybridised 4s4p3d and 5s5p4d basis was used 
in the calculations to describe the Co atoms. Exchange cor- 
relation was taken into account in the framework of the local 
spin density approximation in the form proposed in Ref. 
The lattice constants were a = 2.5071 A and c =4.0695 A. 
38400 k-points were used to sample the Brillouin zone. Re- 
sults are shown in Fig. |2] They are in good agreement with 
previous theoretical results^i. In the output of the calcula- 
tion we find that tensor elements that should be zero due to 
symmetry are of the order of 10 9 s~' while we find a differ- 
ence between a yy and <5 ZZ of the order ofl0 13 ..10 12 s~'. Thus 
the symmetry of our calculated tensor is in agreement with 
Eq. \26\ and Eq. d30i . We conclude that the calculated dif- 
ference between a yy and a zz is a signature of second order 
magnetic contribution. However if we change numerical pa- 



B. Optical response of the single crystal 

We have calculated the optical response in polar MOKE 
geometry with perpendicular incident light with our transfer 
matrix approach. We find that the optical response depends 
strongly on the direction of the polarization vector in the sur- 
face plane. If the polarization vector is along one of the main 
crystal axis, birefringence is absent and the optical response is 
similar to common polar MOKE. When the polarization vec- 
tor is turned away from the main crystal axis the optical re- 
sponse is a combination of crystallographic birefringence and 
a magnetic effect. We find that birefringence starts to be im- 
portant at about 3°. Results are shown in Fig. [3] Directions 
of the polarization vector are in one quarter of the full circle 
in the surface plane which is choosen symmetrically around 
the crystallographic x-axis. For directions of the polarization 
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FIG. 3: Calculated optical response of { 1 120) hep Co in polar MOKE 
geometry. The polarization vector is parallel to the crystallographic 
x-axis at zero angle. Curves shifted to higher values just below 5 eV 
correspond to positive angles. 



vector chosen around the crystallographic y-axis, results are 
identical on the scale of the plot. The latter is a non-trivial 
result. Since the crystallographic x- and y-directions are dif- 
ferent one would expect independent results in half of the full 
circle. It can only be understood by stressing that the birefrin- 
gence is large compared to the magnetic effect (see below and 
Sec. IVID> . The solid curves show the case when the polar- 
ization vector is parallel to a main crystal axis. The optical 
response is similar to the polar MOKE of hep (0001) Co^iSi. 
To a good approximation it can be regarded as a common po- 
lar MOKE response without birefringence. The dashed and 
dotted curves show the optical response for cases when bire- 
fringence is important. If the polarization vector has an angle 
of ±5.5° relative to the main crystal axis the birefringent con- 
tribution has about the same magnitude as the magnetic effect. 
It reaches its maximum at an angle of ±45°. At this angle it is 
about one order of magnitude larger than the magnetic effect. 

The present system has been investigated experimentally in 
detail by Weller et ali. In this experiment different samples 
were used at least one of which was a polycrystal with two 
types of crystallographic domains related to each other by a 
90° rotation around the surface normal. Experimental results 
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FIG. 4: Calculated optical response of bicrystalline hep (1120) Co. 
At 0° polarization vector is along a main crystal axis of one of the 
domains. 



do not report birefringent contributions nor a dependence on 
the direction of the polarization vector. So our theoretical re- 
sults for the single crystal presented here are very different 
from experimental findings. Still, there is no direct disagree- 
ment between theory and experiment simply because it is pos- 
sible that the experimental data that was taken actually corre- 
sponds to the case when the polarization vector is along a main 
crystal axis. For this case there is good agreement with theory 
(see Fig. [5}. However we believe that this is not what was 
happening. Rather we speculate that during measurements 
at some point different directions of the polarization vector 
were used and still basically common polar MOKE was found 
without substantial dependence on the direction of the polar- 
ization vector. Let us for the moment focus on the sample 
which we know is a polycrystal. Then the conclusion is op- 
tical response of a polycrystal with two domain orientations 
is fundamentally different from the optical response of a sin- 
glecrystal, so in order to describe experiment correctly it is 
important to consider the full polycrystal rather than a single- 
crystalline sample. 



C. Optical response of the bicrystal 

We have calculated the optical response of a polycrystal 
with two domain orientations. Our approach was to calculate 
average Stokes parameters from our transfer matrix calcula- 
tion as described in Sec. [II]] Experimental data about the dis- 
tribution of domain sizes and intensities shining on them was 
not known so we had to make an assumption here. We expect 
that crystal growth occurrs with equal preferrence in both of 
the two domain orientations so total surface areas should be 
the same and total intensity of the incident light should be de- 
vided equally among the two orientations. 

Results are shown in Fig. |4] In general, we find now for 
any direction of the polarization vector that our calculated op- 
tical response is similar to common polar MOKE and theoreti- 
cal results are now in good agreement with experimental data. 
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The birefringent contribution, which for the single crystal was 
the dominant contribution to the optical response, is now av- 
eraged out. However birefringent contribution is averaged out 
completely only in the ellipticity (in our computational result 
variation under change of the direction of the polarization vec- 
tor is of the order 10~ 4 °) while in the rotation it is still present. 

In general the results are quite surprising: For the single 
crystal birefringence was about 10 times larger than the mag- 
netic effect. For the polycrystal it is averaged out so strongly 
that it is now smaller than the magnetic contribution. How is 
this possible only due to the presence of one additional do- 
main orientation? And secondly: why is the birefringent con- 
tribution completely missing in the ellipticity but still present 
in the rotation? It is important to find out the general mecha- 
nism behind this. 

We have considered average Stokes parameters for poly- 
crystals with ordered domains analytically. We find that the 
optical response strongly depends on the in plane symme- 
try of the domain orientations. In the majority of cases, or- 
dered polycrystals are equivalent to polycrystals with random 
domain distribution and thus optical response is independent 
of the direction of the polarization vector. In particular we 
can prove that the Stokes parameters Sq and S3 are identical 
to those of a random polycrystal if and only if the in plane 
symmetry of domain orientations is larger than 2-fold and the 
Stokes parameters S\ and S2 are identical to those of a random 
polycrystal if and only if the symmetry of domain orientations 
is not 1, 2 or 4. The prove is given in the appendix. Analytical 
findings are in good agreement with the computational result 
we present here for the hep (1120) polycrystal with two do- 
mains. In particular they explain the different behaviour of 
averaging out in rotation and ellipticity (only S\ and 52 enter 
in the rotation, Eq.^] while mainly S3 enters in the ellipticity, 
Eq-El now note the polycrystal with two domains oriented by 
a 90° rotation has 4-fold symmetry). The analytical findings 
have an important consequence for experiments. They imply 
that if only a few ordered domains are present inside the illu- 
minated area the optical response will always be very close to 
common polar MOKE. 

We now conjecture that the second sample that was investi- 
gated in experiments (the Ru(l 120) sample) was also a poly- 
crystal (the presence of few ordered domains in the illumi- 
nated area is enough). For any direction of the polarization 
vector that was possibly considered in experiment we imme- 
diately have agreement with theory. Summarizing comparison 
of theoretical data with experiment is shown in Fig. [5] Data 
for (0001) hep Co are shown for comparison. The theoretical 
data for (0001) hep Co has been calculated in the same way as 
the data for (1 120) hep. It is in good agreement with previous 
theoretical data 2 ^2i. 



D. Anisotropic polar MOKE 

The goal of the previous experimental work of Weller and 
coworkers was to find a manifestation of magneto-crystalline 
anisotropy in the magneto-optical response. They investigated 
how the optical response changes when the relative orientation 
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FIG. 5: Optical response of hep ( 1 120) Co in polar MOKE geometry. 
Theoretical data is for special cases of the direction of the polariza- 
tion vector. Experimental data is due to Weller et al. Results for hep 
(0001) are shown for comparison. 



between magnetization and crystal lattice is changed while the 
polar measuring geometry as well as other parameters of the 
experiment are kept (lattice parameters, crystal growth qual- 
ity, etc.). It was found that the optical response of hep (0001) 
and hep (1120) is different. These results were explained by 
the dependence of the absorptive part of the refractive index 
on the angle between crystallographic c-axis and spin mo- 
ment. 

We know now that the electrodynamics part of the problem 
is much more complicated. It is important to calculate the full 
optical response including crystallographic birefringence and 
also the polycrystalline nature of the sample has to be taken 
into account. So it is important to check if the main conlucions 
given in the experimental work stil hold. As we will see below, 
the answer is yes. 

From a theoretical point of view the situation is the follow- 
ing: We have common polar MOKE in hep (0001) (no bire- 
fringence, optical response is independ of direction of polar- 
ization vector) and a combination of birefringence and mag- 
netic response with strong averaging out of birefringence in 
the hep (1120) polycrystal. So optical responses are funda- 
mentally different. Nevertheless in both systems the magnetic 
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contribution to the optical response originates from the tensor 
element e™. We would say that we have measured anisotropy 
in the magneto-optical constants if we can conclude from the 
measurement that e^. has changed due to change of the mag- 
netization direction. So what we want to show now is that the 
difference in magneto-optical response between single crys- 
talline hep (0001) and poly crystalline hep (1 120) is basically 
only determined by the change in 8xy. Admittedly we do not 
think this can be proven rigorously, however what we can do 
is to calculate the optical response of the hep (1120) crys- 
tal with a dielectric tensor were we substitute e vy by E n and 
vice versa. We can also use the average i (e XT + e JV ) for both. 
We find in any case the optical response is very close to both 
the result obtained for the single crystal with the polarization 
vector along a main crystal axis and for the polycrystal. All 
these cases are much closer to each other than to the result 
for hep (0001); see also Fig. [5] The conclusion is that the 
difference between optical response of hep (1120) and hep 
(0001) is mainly due to the change in e^. In this sense it may 
be regarded as anisotropic polar MOKE or a manifestation of 
magneto-crystalline anisotropy in the optical response. 



VII. POLAR MOKE OF CR0 2 
A. Optical conductivity 





We have calculated the optical conductivity tensor of 
(010) CrC>2 with the first principles approach as described 
in Sec. II VI The basis set was constructed from 4s4p3d and 
AdAf, (respectively 2s2p and 3s3p) orbitals for the chromium 
(respectively oxygen) sites. The lattice constants and position 
parameters were a = 4.421 A, c = 2.916 A and x = 0.3053 as 
it was used in Refs. Il2l26l27l28l 32768 k-points were used to 
sample the Brillouin zone. Exchange correlation was treated 
in the same way as in the calculation for Co above. 

The magnetic moment per CrC»2 m = 2.0/jb and total en- 
ergy per unit cell as we ll as the DOS agree well with those 
given in Refs. 12 18 28. Fig.|6]shows our calculated optical 
conductivity tensor. Results are in good agreement with pre- 
vious theoretical findings 



B. Optical response of the polycrystal 

If thin films of Q-O2 are deposited on single-crystalline 
AI2O3, polycrystalline growth is observed. Crystallites order 
6-fold symmetrically with an a-axis oriented perpendicular to 
the surface 29 . Experimental results suggest that the sizes of 
crystallites in such films are typically of the order 0.1-10/jm. 
For the lower limit we are in a regime where interference ef- 
fects start to play a role. Consequently the optical response is 
no longer a purely incoherent wave and can in general not be 
described by average Stokes parameters. We exclude this case 
here. For the upper limit the optical response of the polycrys- 
tal is well described by average Stokes parameters. 

We have calculated the optical response of polycrystalline 
(010) Cr02 with 6-fold symmetric domain ordering. We find 



to (eV) 

FIG. 6: Calculated optical conductivity tensor of (010) Cr0 2 . 



that the optical response is independent of the direction of the 
polarization vector. Results are shown in Fig.0 They are in 
good agreement with experimental data. 

Also, results are in good agreement with analytical find- 
ings given in appendix|X] the 6-fold symmetric polycrystal is 
a member of the isotropic class which implies that crystallo- 
graphic birefringence is averaged out completely both in the 
rotation and in the ellipticity. 

Results have an important implication. In a previous theo- 
retical work Uspenskii et al. 12 derived an approximative an- 
alytic expression for the polar MOKE of a polycrystalline 
surface with two-dimensional random domain distribution. It 
reads 



¥+& = 



2e x 



(33) 



Here the roots are taken in the upper complex half plane. 
From the more recent experimental works^it is clear that do- 
main distribution of polycrystalline QO2 is actually not ran- 
dom rather it has 6-fold symmetry. So Eq. d33l is in general 
not applicable. However, from the analytical results of ap- 
pendix|A|we know now that the optical response of the 6-fold 
symmetric polycrystal is equivalent to the optical response of 
a random polycrystal of the same material. Thus, the valid- 
ity of the approximative expression is extended to the whole 
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0.2 - y\ % 




to (eV) 



FIG. 7: Calculted optical response of a-axis textured 6-fold sym- 
metrically ordered CrC>2. The solid line shows the rotation l|/ of the 
polarization ellipse. The dashed line shows its ellipticity %. 

isotropic class. Hence, indeed the optical response of CrC>2 
can be calculated with Eq. J33I . 

We have calculated the optical response also with the ap- 
proximative expression. Results differ from the rigorous result 
obtained with our transfer matrix calculation and subsequent 
determination of exact average Stokes parameters in the fourth 
relevant digit. This shows that (for &O2) the approximative 
expression is actually very good. Also it shows that compu- 
tational results are in very good agreement with the rigorous 
analytic treatment given in the appendix. 



VIII. SUMMARY AND CONCLUSION 

We have calculated the polar magneto-optical Kerr effect 
for hep (1120) Co and for (010) C1O2. Our approach was 
based on first-principle calculations of dielectric tensors. We 
have addressed the electrodynamics part of the problem, i. e., 
the extraction of MOKE from dielectric tensors, with a trans- 
fer matrix method. We could describe simultaneous occur- 
rence of birefringence and magnetic effect that is present in 
the systems. For polycrystals average optical response was 
described by exact average Stokes parameters taking into ac- 
count the real orientations of domains. 

For hep (1120) Co we found that a single crystal optical 
response depends strongly on the direction of the polariza- 
tion vector. If the polarization vector is along one of the main 
crystal axis optical response is very similar to common po- 
lar MOKE and moreover for the two crystal axis the optical 
response is basically the same. If the polarization vector devi- 
ates more than about 3° from one of the main crystal axis bire- 
fringence is important. For larger angles it dominates over the 
actual magnetic effect. To explain experimental data we had to 
stress that samples investigated in experiment were polycrys- 
tals. We could show that already the presence of two domain 
orientations leads to a strong reduction of birefringent con- 
tribution in the magneto-optical response. Finally we could 
show that the previous interpretation of experimental data in 



terms of a manifestation of magneto-crystalline anisotropy in 
the optical response remains valid. 

For poly crystalline (010) Cr02 we found that the birefrin- 
gent contribution to the optical response is averaged out com- 
pletely. We could verify that a previous approximative an- 
alytic expression describes the optical response exactly also 
for the case of realistic domain orientations. 

The results of our LDA calculations for both hep Co and 
Cr02 are in very good agreement with the experimental data 
(assuming that the data for Co are for a bi-crystal). This is not 
trivial since, in general, correlation effects might be essential 
for the electronic structure of transition metal ferromagnets 30 . 
The effect of local Coulomb interactions on magnet o-optic al 
properties of Fe and Ni has been calculated in Refs. l3ll32l in 
a framework of dynamical mean-field theory (LDA+DMFT 
approach). It appeared that, whereas for Ni the correlation 
effects are important, for Fe there are almost no difference be- 
tween LDA and LDA+DMFT results for optical and magneto- 
optical properties. Our results show that probably correlation 
effects are not very important also for magneto-optical proper- 
ties of Co. As for ferromagnetic Q-O2 recent analysis 33 shows 
that it should be considered rather as a weakly correlated sys- 
tem so a success of our calculations is not surprising. 
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APPENDIX A: CLASSIFICATION OF POLYCRYSTALLINE 
SURFACES 

Most polycrystalline surfaces occurring in nature have ei- 
ther a three-dimensional distribution of domain orientations or 
a two-dimensional distribution with only few domain orienta- 
tions that are related to each other by a rotation round the sur- 
face normal. Three-dimensional distribution is found for sur- 
faces of bulk polycrystals such as, e.g., natural iron. Ordered 
two-dimensional distribution is often found when thin poly- 
crystalline films are grown on single-crystalline substrates. In 
the case of three-dimensional distribution, the domain orien- 
tations are often to a good approximation random. The aver- 
age polar MOKE of a three-dimensional random polycrystal 
is obviously independent of the direction of the polarization 
vector in the surface plane. We skip this case here as well 
as other three-dimensionally ordered polycrystals. Rather we 
focus on polycrystals with a two-dimensional distribution of 
domain orientations. We call a surface n-fold symmetrically 
ordered if the crystallographic structures of all domains can 
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be mapped onto each other by an «-fold rotation around the 
surface normal. We will also use a notion of two-dimensional 
continuously distributed polycrystalline surface. By that we 
mean a polycrystalline surface in which the crystallographic 
structures of the domains can be mapped onto each other 
by suitable continuous rotations around the surface normal 
and all possible orientations occur. This corresponds to two- 
dimensional random domain orientations. Also for this case, 
the average polar MOKE is obviously independent of the di- 
rection of the polarization vector. 

We show now that for most polycrystals with symmetrically 
ordered domains the average polar MOKE is equivalent to the 
average polar MOKE of a continuously distributed polycrystal 
of the same material. 

In particular we prove the following statement. The average 
Stokes parameters (So) and (S3) are identical to those of a 
continuous polycrystal if and only if the in plane symmetry 
of domain orientations is larger than two-fold and the Stokes 
parameters (Si) and (S2) are identical to those of a continuous 
polycrystal if and only if the in plane symmetry of domain 
orientations is not 1, 2 or 4. 

We begin the proof by considering the light reflected from 
a single domain. If reflection is described by means of trans- 
fer matrix method, then, for any wave vector and frequency, 
the complex amplitude of the reflected wave is a linear map- 
ping of the complex amplitude of the incident wave. This can 
be seen directly from the main linear equation, Eq. dl4> . Fur- 
ther, in case of normal incidence, the incident and the reflected 
amplitude vectors may be represented in a common coordi- 
nate system parallel to the surface plane. Thus, if E m and 
E re f l are respective 2-vectors, there is a linear transformation 
T : C 2 -> C 2 such that 



E refl = j E h 



(Al) 



Now let some other domain be identical to the previous one 
up to a rotation 



*(q>) = 

around the surface normal. 

?refl 



coscp sincp 
-sincp coscp 



(A2) 



If E 2 is the amplitude of the wave reflected from the sec- 



ond domain, we have 

jprefl 



:R(y)TR-\y)E 



Now consider an n-fold symmetrically ordered polycrystal. 
Introducing angles % = 2tc|, k= 1, ... ,n, the amplitude vec- 
tors 2?I of the reflected waves are 



Ef =R(y k )TR- l (y k )E" 



(A4) 



By Eq. dl5> the average Stokes parameters are 

(S) = £S ^R(2%^)TR- 1 (2%^-)E m ^ . (A5) 
For a continuously distributed polycrystal we have 

2rc 

(S) = i Js(R(<?)TR- l (<?)E'") dep. (A6) 



It is shown in appendix [C] that the latter two expressions 
are equal in the first and last component if and only if 
n { 1 , 2} and in the second and third component if and only 
if n £ {1,2,4}. This finishes the proof. 

The latter statement is fundamental for the understanding 
of the average polar MOKE of polycrystals. It naturally de- 
cides thin polycrystalline films into three classes: 2-fold sym- 
metrically ordered films, 4-fold symmetrically ordered films 
and all others including two-dimensional random orientation. 
Further, it implies that for polycrystalline films out of the first 
two classes the optical response does in general depend on 
the direction of the polarization vector, thus the birefringent 
contribution to the optical response is in general not averaged 
out. On the other hand it implies that for polycrystals out of 
the last class optical response is independent of the direction 
of the polarization vector, thus the birefringent contribution to 
the optical response is averaged out. 



APPENDIX B: SYMMETRIC SUMS OVER POWERS OF 
TRIGONOMETRIC FUNCTIONS 



We prove a statement about symmetric sums over powers 
of cos and sin. 

Let / : R -> K and q e Q. Then the identity 



271 



(A3) holds for pairs /, q 



cos 



sin 



1 

2 

-2 I 

' 2 

cos sin , 



4 , 2-2 

cos + cos sin 



• 4 , 2-2 

sin +cos sin 



13 1 

2 '4 + 8 
1 3 1 

2 '4 + 8 



(Bl) 



cos~ sin + shr cos , 

• 4 4 rv 

sin — cos , 
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if and only if n { 1 , 2}, and for pairs /, q 

4 1 3 
cos , - • — 
2 4 

.4 13 
sin , - ■ - 

' 2 4 



2 • 2 

cos sin 



cos~ sin , 
sin 3 cos , 



1 



(B2) 



2 ■ 2 

cos sin 



■ 4 2-2 

sin — cos sin 



3 • -3 
cos sin — sin cos 



1 3 

2 "A 
1 3 
2~'4" 





if and only if n £ {1,2,4}. 

We begin the proof by considering sums of the form 



1 



m g N, 



(B3) 



where n e N, n > 2. 

If m = ln with some I e N, then 

in in 



1, 



whereas if n = I m with some leNwe have 



1 



„mi2ni 



I in 



e 



„<2tc4 



1 



/ 



4=1 



— V e'^' 1 ' = — m 



„/2ji 



7 =0. 



(B5) 



Now let /? and g be the largest prime numbers occurring 
in the prime factorisations of n and m respectively. Let 
F ; , = { {0 , 1 , . . . , p — 1 } , • , +} be the prime field of the mod- 
ulo classes of p in the common sense. Let m ■ F C N be the 
set {0 , lq\p , 2q\p , . . . , (p — 1 )q\p}- We devide the set 
of complex numbers occurring in Eq. ( IB3I into s subsets. We 
chose s such that n = s- p and consider 



> = { 



A,= 
Then 



i2n^+mi2%i 



,k€m-¥ p 



z€At 



(B6) 



E 



<i2xi 



kem-¥„ 



= if q < p. 

(B7) 



Using Eq. (IB4> . Eq. (IB5> and Eq. (IB7I . we can calculate the 
sum given by Eq. JB3I with some m G N for any n £ N. We 
consider the cases m = 2 and m = 4. 

For m-2we obviously have a largest prime factor q = 2, 
i.e., by Eq. (IB7> the sum vanishes for 

n = 3, 6,7,9,10,11, 12,... 



and any other natural number containing a prime greater than 
or equal to three in its factorisation. For n = 4, 8, 16, . . ., the 
sum vanishes by Eq. ( IB5> . while for n = 1 ,2 the sum is one 
by Eq. ( IB4I . Thus we have 



-f e 2i2, »=0 if and only if ng{l,2}. 
"i i 



(B8) 



For m = 4, we have a largest prime factor g = 2 as well, i.e., 
the sum vanishes again for 

« = 3,6,7,9, 10, 11, 12, . . . 

and any other natural number containing a prime greater than 
or equal to three in its factorisation. For n = 8, 16,32, . . ., the 
sum vanishes by Eq. (IB5> . while for n = 1,2,4, we obtain by 
Eq. (IB4> that the sum is one. Thus 



1 " k 

-Ye 4 ' 2x T,=0 if and only if {1,2,4}. (B9) 

n k=l 

We prove the first line of Eq. JB It . The identity 
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-Ly c os 2 ((p)<i(p=i (B10) 
o 

(B4) follows from the more general formula 



it/2 



sin 2a+1 ((p) cos 2 P+ »(c P )^cp= r (° t + 1 ) r (P + 1 ) (BU) 
m VHV Y 2r(a+l + P + l)' 



where T is the gamma function 34 . Further, 

fa** +e -2a*t 

n 4^L 

1 



n 111" 
^cos 2 (2^) = i + i.ii: 



(B12) 



if and only if n { 1 , 2} , 



where we have used Eq. (IB8> for the last line. 

In a similar way, the second and third line of Eq. JBH fol- 
low from Eq. ( IbTTI and Eq. dB8l . 

Next we prove the first line of Eq. ( IB2I . Once again, we 
refer to Eq. JB1U to see that 



2n 

— / cos 4 ((p)J(f>= 



27C 



2 4 



(B13) 



On the other hand 
I£cos 4 (2^) 

" k=\ 

13 11 

~ 2 '4 + «' 16 



„4/27[i 



-4e 



-4e 



— 2i2nS 



1 3 

-•- if and only if B {1,2,4}, 



(B14) 
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where we have used Eq. (IB8> and Eq. ( IB9> for the last line. 
In a similar way, we find 



1 " 1 t, 

- Vsin 4 (27^) = ---, 



i£cos 2 (2^)sin 2 (2^) = i 
"fc=i 



(B15) 



(B16) 



which expresses, that in a total intensity measurement, the re- 
flected light of a polycrystal illuminated with a single incident 
beam is not distinguishable from the reflected light of a sin- 
gle crystal illuminated with several beams with respective ori- 
entations of the polarization vectors. We denote c = cos((p), 
s = sin(cp) and 



R 



c s 
—S c 



(C4) 



and 



1 " 

- £cos 3 (27i;£)sin(27i;£) = (B17) 



k=l 



1 " 

- £cos(27i;£)sin 3 (27i;£) = (B18) 

n k=i 

if and only if n ^ {1,2,4}. This gives the first five identities 
ofEq. dH3. 

To see that the last four lines of Eq. (IB It hold, we add the 
corresponding expressions obtained above and find that in all 
cases the sums over fourth powers cancel, while sums over 
second powers remain. In contrast to that, also the fourth or- 
der sums remain in the expressions for the last three lines of 
Eq. (IB2t . This finishes the proof. 



Thus 

tr 1 — I 

y CTy X +STyy STy X +CTyy 

Introducing £ = (a,b), we have 



TR l _ l cT xx + sT xx . ~sT xx + cT r 



(C5) 



[TR-'E] x [TR- l E] x = [(cT^ + sT^a + i-sT^ + cT^b] 



[{cT xx + sT xy )a + (sT xx + cT^b] 



c 2 T xx T xx + csT^Tjcy + csT„T. 



T I 2 7 1 t 1 



+ (-csTxxTja + c^-T^T^. - s 2 T xy T xx + csT^T^) . 
+ (scTxxTxx - s^xxTxy + (P'TxyTxx + csT^T. 
~\~ (s T XX T XX scT xx T X y scTxyT xx -\-c T xy >T x 



(C6) 



APPENDIX C: AVERAGE STOKES PARAMETERS FOR 
n-FOLD ROTATED 2x2 LINEAR TRANSFORMATIONS 

LetZs e R 2 , R(q>) : K 2 — > K 2 be a rotation by an angle (p and 
T : C 2 — * C 2 a linear transformation of most general symme- 
try. Let Sj : C 2 — > M, j = 0, 1,2,3 be the Stokes parameters. 
Then 



l tSj(H2^,)TR(2K^- l E) 



■ k=l 



2% 



(CI) 



= —J Sj(R(<p)TR-\ 9 )E) dip 



holds for j = 0,3 if and only if n ^ {1,2} and for j = 1,2 if 
and only if n g {1,2,4}. 

We treat the four cases separately. 
Consider So- 

First we evaluate the expression for the Stokes parameter 
occurring in Eq. ( ICU . Dropping the angular argument of the 
rotation, we get from Eq. dl5> 



S Q (RTR- l E) = [RTR- [ E] x [RTR- l E] x 



+ [RTR- l E] y [RTR- l E] y (C2) 
= (RTR~ l E , RTR~ l E), 

where (•,•) denotes the standard scalar product in C 2 . R is 
orthogonal, thus 

(RTR~ l E , RTR~ l E) = (TR~ l E , TR~ l E), (C3) 



Thus by the first three lines of Eq. IB II 

2j[ 



1 /[77? '((p)/. i ■'< 'f«p)£] i( /«p 



(C7) 



1 # 



-Y J [TR- l (2K^E} x [TR-^2^)E] 



k=l 



if and only if n ^ {1,2}. From Eq. IC5I we see, that the y- 
component of the transformed vector has the same form as 
the jc-component. This completes the proof for So. 

Consider Si and $2. 

In contrast to So, both Si and S2 are no scalar products. 
Thus, we have to evaluate the full expression RTR~ l E. With 
Eq. (|C4) and Eq. C3 we get 



RTR 1 = I cljxx + CsTxy + CsTyx + s2Tyy 

\ -CST XX - S 2 J„ + C 2 T yX + CSTyy 



-csT xx + c 2 T xy - s 2 T yx + csT yy 
s 2 T xx - csT xy - scT yx + c 2 T yy 



(C8) 



We denoted = (a,b) as before andRTR l E =E' = (E' x ,E' y ) 
Then 



E' x = (c 2 T xx + csT„ + csT yx + s 2 T yy )a 

+ {-csT xx + c 2 T xy - s 2 T yx + csT yy )b, 

E' y = {-csT xx - s 2 T xy + c 2 T yx + csT yy )a 

+ {s 2 T xx - csT^ - csT yx + c 2 T yy )b. 



(C9) 
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By Eq. d!5l >. we have With Eq. d!5l > we obtain 



S x =E' x E' x -E' y E' y (CIO) 



and 



S 2 =E' x E' y + E' x E' (Cll) 



If we evaluate these expressions by substituting Eq. dC9> . ev- 
ery resulting term contains factors of cos and sin of the form 
considered in Eq. JBH or Eq. JB2I >. These terms might add up 
to combined terms out of the last four lines of Eq. iBli or the 
last three lines of Eq. ( IB2t . To check, if we have at least one 
independent term out of Eq. JB2ll . we focus on expressions 
with a factor T XX T XX . We obtain 

S\ =T XX T XX (c 4 a 2 - c 3 sab + c 2 s 2 b 2 - c 2 s 2 a + cs 3 ab - s 4 b 2 ) 



S-i=i(E'E'-E' x E') 
--T xx T xx ({c 4 -c 2 s 2 )a 2 ~(s 2 c 2 -s 4 )b 2 + (cs 3 -c 3 s)ab) ' 

+ ... (C12) 



[ (c 4 + c 2 s 2 )T xx T xx - (c 4 + c 2 s 2 )T yx T xx 
+ (c 3 s + cs 3 )T xx T yy - (c 3 s + cs 3 )T yy T xx 

and - (/ + c V ) T yy T„ + (s 4 + c 2 s 2 ) T xy T yy 

_ - 

Si =2T XX T XX (c 3 sa 2 + 2c 2 s 2 ab - cs 3 b 2 ) + .... (CI 3) 



+ (c 3 s + rv ; :/ n y rv - (c 3 s + C5 3 )r„r„] a- 

+ [ 2[{c 3 s + cs 3 )T yx T xx - (c 3 s + cs 3 )T xx T y} ] 

Thus, we have independent terms out of Eq. dB2l which can- + 2{(c 3 s + a 3 )T T — (c 3 s + cs 3 )T T ] 

not be combined to terms out of Eq. iBl\ . This proves the _ _ 

cases 5 1 and 52. + ( c JrC s )TxxT yy —(c s +c )T yy T \ x 

Last consider 5 3 . + ^ + " f* + 

Using the same notation as before and the results obtained + (c — s )7iy7^ ;c + (s 4 -c 4 )T yx T xy \ab 

inEq.Swehave +[ (cV+^T,,- (s 4 + c 2 s 2 )T yx T xx 

+ {cs 3 + c 3 s)T yx T xy - {c 3 s + cs 3 )T xy T yx 
+ (cs 3 + c 3 s)T xx T yy - (c 3 s + cs 3 )T xx T yy 
+ (c 4 + c 2 s 2 )T xy T yy - (c 2 s 2 + c 4 )T yy T yx ] b 2 . 

a 2 



= [ c sTxxTxx~ 


C S T X x T xy 


+ c 4 T xx T yx + 


C STxxTyy 


C S~TxyT XX 


Txy T xy 


+ c 3 sT xy T yx - 


~ C S Txy Tyy 


C S Ty X T XX 


-cs 3 T yx Txy 


+ c 3 sT yx T yx - 


~C S Ty X Tyy 




-cs 3 T yy T xy 


+ c 3 sT yy T yx - 


~ C S TyyTyy 


H~ [ c s T X xT xx 


c sT xx T xy 


-c 3 sT xx T yx - 


~ C TxxTyy 


~\~CS TxyT XX ~ 


- C S T X y T xy 


-c 2 s 2 T xy T yx 


+ C STxyTyy 


-\-cs Ty X T XX ~ 


- C S Ty X T X y 


C Ty x Tyx 


+ C STy X Tyy 



+S 4 T yy T xx ~ rV ; '/;,//\,, - CS 3 TyyTy X + C^SpTyyTyy 

~r~c s T XX T xx -+- c sTxxT xy c sTxxTy X — c s T xx Tyy 

-C 3 STxyTxx - C 2 S 2 TxyT X y + C^TjcyTyx + c 3 sTxyT yy 

~\-CS Ty X Txx~\~$ TyxTxy C S Ty X Tyx CS Ty X Tyy 

-c 2 s 2 T yy T xx - cs 3 T yy T xy + c 3 sT yy T yx + c 2 s 2 T yy T yy ] ab 
- cs T xx T xx -\~c s T XX T ' xy + c s T xx T yx c sT xx T yy 
+c s Tj-yZ 'xx c sTxyT 'xy c sTxyT yx -\-c T^T yy 
s T yx T xx -\rCs T yx T xy ~\~ cs T yx T yx c s T yx T yy 



cs 3 T yy T xx c s TyyTxy c s T yy T yx + c sT yy T yy ]b . j ndeed a jj terms conta j n a factor of cos and sin out of those 

given in Eq. iBl\ . This finishes the proof. 
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